EPJ manuscript No. 

(will be inserted by the editor) 



Adiabatic quantum pumping in graphene with magnetic barriers 

E. Grichuk and E. Manykin 

National Research Center "Kurchatov Institute", Kurchatov Sq. 1, Moscow 123182, Russia 
the date of receipt and acceptance should be inserted later 

Abstract. We study an adiabatic quantum pump effect in a two terminal graphene device with two os- 
cillating square electric barriers and a stationary magnetic barrier using the scattering matrix approach. 
The model employs the low-energy Dirac approximation and incorporates the possible existence of a finite 
band gap in graphene spectrum. We show that in this case valley-polarized and pure valley currents can 
be pumped due to the valley symmetry breaking. For a 5-function magnetic barrier we present analytical 
expressions for bilinear total and valley pumping responses. These results are compared to numerical ones 
for a double J-function, a square and a triple square magnetic barriers. 
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1 Introduction 

Graphene, a single atomic layer of carbon atoms packed 
into a honeycomb lattice, is now being intensively studied 
both theoretically and experimentally. In an endeavor to 
construct new nanoelectronic devices based on graphene, a 
particular attention of researchers is paid to its unconven- 
tional and somewhat counterintuitive electronic proper- 
ties. Graphene electronic spectrum comprises two valleys, 
referred to as K and K' , and in the low energy approx- 
imation electron dynamics in each valley is governed by 
a 2D Dirac equation for massless particles. Thus, quasi- 
particles in graphene have linear spectrum with the Fermi 
velocity vp ss c/300. This results in a peculiar behavior 
of electrons in this material, which manifest itself in the 
Klein tunnelling [H[5], the half- integer quantum Hall ef- 
fect [3], the Veselago lensing of electrons [4] and many 
other effects 0IMZli]. 

When graphene is considered as a basis for nanoelec- 
tronics (e. g., graphene FETs), the Klein tunneling through 
a potential barrier is of special interest. Touching of con- 
duction and valence bands and the existence of propagat- 
ing modes below the barrier allow electrons to penetrate 
it with high probability (with unit probability at normal 
incidence). Thus, this effect hampers the confinement of 
electrons in graphene-based field effect devices. A num- 
ber of ways have been proposed to circumvent this ob- 
stacle. For example, a graphene sheet can be cut to form 
nanoribbons, nanoislands, nanorings and other nanostruc- 
tures. Operating principles of many suggested graphene 
devices based on such structures crucially depend on the 
atomically precise edge configuration, e. g., a zigzag or an 
armchair edge of graphene nanoribbons. Their fabrication 
is an experimentally challenging task and forces one to 
look for alternative methods. 



Let us briefly mention two of them which are rele- 
vant for the present paper. First, electron confinement 
in graphene can be achieved by means of inhomogeneous 
magnetic field [ WTUIITT] that can be generated, in par- 
ticular, by depositing ferromagnetic strips (gates) on the 
top of a dielectric layer covering a graphene flake. Sim- 
ilar technique is now widely utilized to create magnetic 
barriers in structures with conventional two-dimensional 
electron gas. Their graphene counterparts are expected to 
be realizable in the near future with various useful appli- 
cations, for example, in graphene-based spintronics. Sec- 
ond, a special choice of a substrate can break graphene 
sublattice symmetry by introducing staggered sublattice 
potential and induce a finite energy band gap [T2j[T3] . De- 
pending on a substrate (e.g., ft,-BN or SiC), values from 
a few tens of meV up to a few tenths of eV have been 
reported. The combination of an external magnetic field 
and an induced band gap can be used to generate and 
detect valley-polarized currents in graphene — a necessary 
step toward the so-called graphene "valleytronics" (valley- 
based electronics) , which exploits peculiar distributions of 
carriers in the valley space [mUMIM!?] ■ 

Most research of electronic transport in graphene is 
focused on stationary problems. A variety of new effects 
emerges when one considers non-stationary ones. An inter- 
esting phenomenon, initially due to D. J.Thouless [TS], is 
a quantum pump effect, in which a periodic modulation of 
parameters of a quantum system produces a finite dc cur- 
rent through it even in the absence of an external bias. 
Quantum pumping in graphene has recently attracted in- 
creasing attention of researchers [TW2ll2l1l2"2"II23[|24l,l251 
26 .27,28,29,30,31, 32 . The unusual electronic spectrum 
of graphene was demonstrated to have a significant im- 
pact on the effect. In particular, the important role of the 
Klein tunneling and evanescent modes was stressed |19[ 
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To model the electron transport in the graphene de- 
vice, we employ the low-energy Dirac approximation. The 
electronic transport is assumed to be completely phase- 
coherent between the leads. The single-valley Hamiltonian 
of the device reads 
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1. (a) A schematic structure of a proposed graphene de- 
vice. The device is formed by a wide graphene ribbon with 
two square electric barriers (produced by top metallic gates) 
of heights Ui and U r and width L as well as (b) a single 6- 
function magnetic barrier or (c) a square magnetic barrier of 
width 2L or (d) a double 5-function magnetic barrier or (e) a 
triple square barrier of width 3L. 



20, 2l]. The potential use of a quantum pump effect in 
graphene-based spintronics was also discussed [23,24,25, 
125112711281125]. 

The interaction of graphene with laser field falls under 
the same category of time-dependent phenomena, and has 
been intensively studied [34,33,35,36,37 . It was demon- 
strated that external radiation field can qualitatively change 



where vp is the Fermi velocity in graphene, p is the canon- 
ical momentum operator, A(x) is the vector potential cor- 
responding to a magnetic barrier and U(x) is the scalar 
potential corresponding to electric barriers. To simplify 
the notation, we put Ti = e = 1. 

For electric barriers we adopt the following profile, 
which is translationally invariant in the y-direction: 



U h -L < x < 0, 
U(x) = \u r , < x < L, 
0, otherwise. 



(2) 



the electronic transport properties of graphene. In particu- 
lar, it can be employed for tuning a band gap in graphene [33] spending vector potential takes the form A(x) 
for generating valley-polarized current in bilayer graphene [33 ,Af ith 



The sharp edge approximation, which is suitable for 
analytical calculations, to physically relevant smooth bar- 
riers is justified if the edge smearing length of barriers is 
smaller than the Fermi wavelength of electrons, being at 
the same time larger than the graphene lattice constant 
to suppress the intervalley scattering. 

To obtain an analytical solution, wc first consider a sin- 
gle 5-function magnetic barrier B(x) = 2BL5(x)z, trans- 
lationally invariant in the y-direction with magnetic field 
perpendicular to the graphene sheet and localized between 
electric barriers (Fig. [Tja). In the Landau gauge the corre- 

A(x)y 



and for analysing transport in graphene superlattices us- 
ing a spatial-temporal duality between static spatially pe- 
riodic electric and magnetic fields and time-periodic laser 
field [35 . It is also interesting to note that laser field can 
induce a topologically nontrivial band gap in (otherwise 
gapless) graphene spectrum leading to the formation of a 
topological insulator state in graphene [37] . 

Motivated by possible applications of a quantum pump 
effect in graphene valleytronics, in this paper, we extend 
previous studies of quantum pumping in graphene by tak- 
ing into consideration the valley degree of freedom of elec- 
trons. 



2 Delta-function magnetic barrier 

The system that we examine is a standard two terminal 
quantum pump device that is formed by a wide graphene 
strip of width W in the (x, y)-plane with two electric barri- 
ers whose heights Ui and U r can be periodically modulated 
in time and one stationary magnetic barrier (Fig. [T]) . Ex- 
perimentally such barriers can be realized by depositing 
ferromagnetic strips on top of a graphene ribbon [38] . In 
this paper we limit ourselves to considering simplest pro- 
files of a magnetic barrier: 5-function barriers and square 
barriers. 



A(x) 




x < 0, 
x > 0. 



(3) 



For later use, we define the following length and en- 
ergy scales inherent to the problem: £g — 1/VB and 
El = vf/L. For a typical magnetic field Bq = 0.1 T 
and Fermi velocity vf = 0.54 eV-nm, we have Lq = 
£b = 81.1 nm and Ei — 6.6 meV. The Zeeman split- 
ting Ez = g^isBo = 1.8 • 10~ 3 El is small and will be 
neglected here (the spin degeneracy factor of 2 is omitted 
throughout the paper). 

The eigenfunctions of the Hamiltonian ([!]) are well 
known 0[6]. Due to the translational invariance of the 
system in the y-direction, in the region a where U (x) = 
U a and A(x) = A a , they can be written as ip a (x,y) — 
e iqy ijj a (x) with 



1p a (x) 



1 



Ak a x 



i 



i 



-ik a x 



(4) 



where r\ = vp(k a + iq a )/(E — U a ). The eigenenergy E is 
given by E = U a ± VF^Jk%, + q^. For electron- like states 
(E > U a ) the first and the second states in Q corre- 
spond to left- and right-moving carriers, respectively. Be- 
low only electron-like states in the leads (E > 0) are con- 
sidered. Note that for finite magnetic field, fej ^ k r , and 
the factor 1 / \fk~a, is included to ensure the unitarity of 
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the scattering matrix. The transverse canonical momen- 
tum q, which is conserved in the scattering process, and 
kinetic momentum q a are related by q a = q + A a . In the 
left (x < 0) and in the right (x > 0) regions we have, 
respectively, 



qi = q — BL and q r = q + BL. 



(5) 



Transverse modes indexed by q are not mixed by the 
scattering, and the scattering problem is solved indepen- 
dently for each mode q. This can be done conveniently us- 
ing a transfer matrix method [39) . In this method potential 
profiles are approximated by piecewise constant functions. 
A transfer matrix of the whole system is then calculated 
as a product of (known) transfer matrices describing the 
propagation of an electron through regions with uniform 
potentials and potential steps. A scattering matrix S can 
be extracted from a transfer one using well-known rela- 
tions. 

The device considered in our study is up-down sym- 
metric, i. e. it possesses a symmetry axis parallel to the 
current direction (the x-direction) . Then, the S'-matrix is 
always symmetric |40j and can be cast in the form 



S = 



r t' 
t r' 



(6) 



where r and t (r' and t' = t) arc reflection and transmis- 
sion amplitudes for electrons incident from the left (right) 
lead. 

Suppose now that electric barrier heights £/; and U r are 
periodically varied in time so that a point (Ui(t),U r (t)) 
traverses a closed contour Cjj in the (Ui, U r ) plane, say, 



Ui(t) = U w - 

U r (t) = U r0 



SUi sin(wt), 
- SU r sin(cLi£ - 



(7) 



with some fixed phase shift <j>. This will lead to a finite 
charge being pumped through a device per each cycle. If 
the pumping frequency is small enough (adiabatic regime), 
the transported charge in the mode q can be expressed via 
a "frozen" scattering matrix S(t) = S(Ui(t), U r (i)), which 
depends on the time t as a parameter [4"Tll42| : 



Q(E,q)-- 
where 

n{E, q ,u h u r 



duidu r II(E, q, Ui, U r ), 



Ua 

Er 



(8) 



^Im 



dr* dr dt* dt 



dUi dU r dUi dU r 
fd^dT_ _ dT d V 
2ir \dUi dU r " WidU r 



El 



(9) 



The expression ([8]) can be significantly simplified pro- 
vided the amplitudes SU a are small (6U a <€. El) so that 
II(E, q, Ui, U r ) is approximately constant within Cjj- Then, 
the response of a quantum pump is bilinear in the ampli- 
tudes SUi and SU r : 



Q(E,q)=n(E,q,Uio,U r0 )A u , A v 



Au 



(10) 



where Au = ir5Ui8U r sin0 is the area enclosed by the 
contour C\j on the (Ui,U r ) plane. 

To characterize the pumping, we numerically calcu- 
late the total pumped charge Q(E) and the pumping re- 
sponse x(E) averaged over the transverse modes 



Q(E) 



Au 
E\ 



J2n(E,q), X (E) = J- £7705,,?), (II) 



where summation extends over N m transverse modes that 
are propagating in both leads at the energy E. For wide 
ribbons the sums can be approximated by integrals. 

A succinct analytical solution is obtained for a device 
with Uio = U r o = 0. For the transmission T = \t\ 2 through 
the (^-function magnetic barrier ^ in the absence of elec- 
tric barriers one gets [TU] 



T(E,q) = 



Ahk. 



(ki + k r 



{2L/e B r 



(12) 



where k r 



•jp — q^, a = I, r. The magnetic barrier 
has a finite transparency only if the energy E exceeds a 
minimum value [S] 



E„ 



E, 



(13) 



otherwise propagating modes do not exist in both leads. 

Taking the derivatives of the scattering matrix S in ^ , 
we end up with the following expression for the pumping 
response: 

n(E,q) = 2 T } qi ^n rL l sin(fc ; L) sin(fc r L)x 



7T (kiL) 2 {k r L) 2 
x sin[(fc/ + k r )L] 



LE 
7T \t^E~l 



T 2 sgn(B) 



sin 2 (k r L) — ~ sin 2 (fc;L) 



(14) 



where sgn(z) is a sign function and a factor of 2 was in- 
cluded to account for the double valley degeneracy (see 
next section). 

The second term in this expression contains the factor 
1/1% = B and gives no contribution in the limit of vanish- 
ing magnetic field. The expression then simplifies to the 
result obtained by E. Prada et al. [5D] 



n(E,q) 



2 2(qL) 2 sm 3 {kL)cos(kL) 
7T (fcL) 4 



(15) 



with k = ki = k r and q = qi = q r . It is readily seen that 
the response vanishes for the modes with q = 0, which 
travel normally to the barriers. These modes are insensi- 
tive to the electric barriers Ui and U r due to the Klein 
paradox and therefore cannot be pumped |20j . This can 
also be seen from ([9]). In the absence of magnetic field 
dT/dUi — dT/dU r , and the pumped current vanishes if 
the transmission T is independent of Ui and U r . 
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The Klein tunneling also manifests itself for a magnetic 
barrier of finite height. The first term in (14) contains the 
product qiq r of the transverse kinetic momenta. Hence, 
this term drops out if the mode propagates in the left or 
in the right lead normally to the barrier. In either case, 
the transmission is insensitive to the height of the corre- 
sponding barrier. The contribution of the second term is 
always finite and indicates the phase coherent nature of a 
quantum pump effect: although one of the barriers does 
not affect the transmission probability of an electron, it 
influences its phase 0T]. If, say, qi = 0, then dT/dUi = 0, 
but dT/dU r remains finite, and for the pumping response 
we have ^ 

n(E, q ) = -^§-, d6) 

where a simple relation dtp/dUi = —L/vp = —I/El f° r 
the mode with qi — was used. 

The spatial symmetry of a device can impose certain 
symmetries on the response function II(E,q). It can be 
easily demonstrated that the Hamiltonian (fTl) satisfies the 
symplectic symmetry SH(B)S~ 1 = H{—B) (operation of 
time-reversal in a single valley), where S — ia y C with 
C the operator of complex conjugation. This symmetry 
implies that @3IHi5] 



59 4 

«4 2 



5g 4 



n(q,B) = n(-q,-B) (17) 

and, therefore, 

Q(B) = Q(-B), X (B)=X(-B). (18) 

These relations are valid for arbitrary Uiq and U r o- If the 
electric barrier is left-right symmetric (Uio = U r o) and 
the vector potential is left-right antisymmetric (A(x) = 
—A(—x)), then the pumping response is also invariant 
with respect to the inversion of q alone: 



n(q,B)=n(-q,B). 



(19) 



This equation is a consequence of the symmetry OH(B)0~ 1 
H(B, U r <-> Ui), where O = <r z R x R y with R x (R y ) the re- 
flection about the x-axis, y — > — y (y-axis, x — > —x) and 
U r O Ui denotes the interchange of the electric barriers 
heights. 

The dependences of the pumping response L7(E, q) and 
the transmission T(E, q) on the canonical momentum q 
and the Fermi energy E are shown in Fig. [2] One can ob- 
serve that the perfect Klein tunneling in the vicinity of 
q = is destroyed by the external magnetic field. As it in- 
creases, an interference pattern becomes visible (Fig. [2ji). 
A similar pattern emerges in the transmission through the 
device with finite Uiq and U r o (not shown). 

In Fig.|3]we plot the averaged pumping response x(-E') 
and the total pumped charge Q(E) for different values 
of the magnetic field. In zero magnetic field the averaged 
response saturates at the value of 0.5 [50]. The magnetic 
field destroys the monotonic behavior and diminishes the 
effect at high Fermi energy. When B is swept from zero 
to B , the response is significantly reduced. At stronger 
magnetic field a similar but weaker tendency is observed. 
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Fig. 2. (a)-(d) Contour plot of the pumping response II(E, q) 
and the transmission T(E, q) for the 5-function magnetic bar- 
rier (|3| in graphene with gapless spectrum as a function of the 
canonical momentum q and the Fermi energy E. In panel (c) 
the contributions of the first (n^) and the second (II^) 



terms in ( 14 1 are shown separately. Dashed vertical lines corre- 
spond to qi — and q r = 0. All distributions are even functions 
of q, and only half of each distribution is shown. Width of elec- 
tric barriers is L — Lq. 




(E-EJE L 

Fig. 3. (a) The pumping response x(E) averaged over trans- 
verse modes and (b) the pumped charge Q(E) for a wide ribbon 
as a function of the Fermi energy E for different magnetic field 
strengths. 



In all curves in Fig. [3] the pumped charge Q{E) is a 
positive function of energy, i. e. the pumping is directed. 
However, by tuning the electric barrier heights Uiq and U r o 
a regime, in which Q(E) becomes a sign-changing func- 
tion, can be achieved. Then, the pumping direction can 
be reversed by varying the Fermi energy — a generic fea- 
ture inherent to quantum pumps. It can be used to gener- 
ate pure spin ^5,26,27,28,29^ and pure valley (see below) 
currents. 



3 Valley-polarized current 

In the above analysis a valley degree of freedom was of 
no importance and was accounted of by a factor of 2. The 
symmetry between the K and K 1 valleys is destroyed in 
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gapped graphene with broken inversion symmetry [14,15, 
116) , and this opens the possibility to pump valley-polarized 
currents. 

The graphene Hamiltonian that describes the states in 
both valleys can be written in various unitarily equivalent 
representations. Choosing a basis (ipA, i>B ) T in the K val- 
ley and (— ip' B ,ip' A ) T the K' valley the Hamiltonian takes 
the so-called valley isotropic form 



77 = v F a[k + A{x)] + Ar z a z + U(x) 



with A the site energy difference between two sublattices 
(determined by the substrate) and t z = r the valley index: 
r = 1 (-1) for the K (K 1 ) valley. 

To obtain an analytical solution we again start with a 
8- function magnetic barrier |3| . Eigenfunctions of Hamil- 
tonian ([20| with uniform scalar and vector potentials are 
given by B with rj = v F{k a +iq a ) / (E — U +At). From an 
eigenenergy E = U a ± yj A 2 + v 2 F (k 2 l + q 2 ) it is seen that 
the spectrum has an energy gap 2A. The minimum en- 
ergy is now E m = E LV / (A/E L ) 2 + {L/l B ) A . Using trans- 
fer matrix method and performing calculations similar to 
those in the previous section, we can find the transmis- 
sion T T (E,q) and the pumping response II T {E,q). 



The expression ( 12 ) for the transmission remains valid 
with ki and k r given by k a = \J(E 2 — A 2 )jv 2 F — q 2 . The 
transmission is independent of the valley index r, and no 
valley-polarized current is produced by applying a bias 
voltage. However, as was discussed in Refs. [I51IT5] . the 
transmission acquires a dependence on the valley index r 
in the presence of a non-uniform electric field. Hence, in 
contrast to the transmission, the pumping response for 
a device with Uio — U r o = is expected to be valley- 
dependent. After tedious but straightforward calculations 
we arrive at 



n r(E,q) = ~ J^n^l Ulr sin{hL) sin(AvL) x 



7T {kiL) 2 {k r Lf 



x sin[(fc; + k r )L — reft] + 



1 / LE 




7T \l B E L 



T 2 sgn(73) 
{hL) 2 {k r Lf 



-j r sin(fc r L) sm(k r L + Ttfi r ) 

r 

7; sin(fciL) sm(kiL - T(f>i 



where 



and 



sin < 



EATsga(B)(ki + k r )L 
E l lilrkik r qiq r L 2 l 2 B 



7a = 1/1 



VFQa 



tanf/)„ = 



k a A 
q a E' 



(21) 



(22) 



(23) 



For vanishing band gap the expression reduces to the 
valley-independent result ( 14 ) of the previous section. In 



the absence of a magnetic field we get the following gen- 



eralization of ([15]): 

2 2 [(qL) 2 + (A/E L ) 2 ] sin 3 (fcL) cos(kL) 



n(E,q) = 



(kiy 



m) 

The presence of a band gap breaks the perfect Klein 
tunneling, and the contribution of modes with normal in- 
cidence (having q = 0) to the pumping response becomes 
non-zero. Integrating (24) at small k 7 for the averaged re- 



(20) sponse x(^) we obtain 



X(E) 



E 2 



2E L VE 2 - A 2 



(25) 



Hence, in gapped graphene the averaged response x(E) 
diverges at small Fermi energy E — > A and the pumped 
charge Q(E) ~ (E 2 + A 2 )/E\ tends to a constant value 
(Fig.§. 

To characterize the currents carried by electrons in dif- 
ferent valleys, we introduce the total II c (E,q) and the 
valley 77 v (E, q) responses 



77,(75, q) = n K (E, q) + n K , {E, q), 
n v (E, q) = n K {E, q) - II K ,{E, q) 



(26) 



and, analogously, the total Q C (E) and the valley Q V (E) 
pumped charges. 

From (21 ) for the valley response we get 



2 / L \ EA T 2 sgn(S) 



(k[ + k r )L 



cos[(fc/ + k r )L] sm(kiL) sin(k r L)- 



{kiL)(k r L) 

+ sin[(fc; + k r )L] cos[(fc r - ki)L] \. (27) 



This expression contains the factor BA, and as ex- 
pected, no valley-polarized current is pumped in zero mag- 
netic field or when the band gap closes. 

From the above equations it is clear that the pumped 
charge Q c is invariant with respect to the inversion of mag- 
netic field, whereas the valley pumped charge Q v changes 
sign: 



Qc(B) = Qc{~B), Q V (B) = -Q v (-B). 



(28) 



These relations are dictated by the time reversal sym- 
metry. The mass term Ar z a z in the Hamiltonian (20) 



breaks the symplectic symmetry 5, but the orthogonal 
time-reversal symmetry TH(B)T^ 1 = H(—B) with T = 
—TyCTyC is preserved (here the Hamiltonian is viewed as a 
4x4 matrix acting in the valley and the sublattice spaces, 
and the operator r y interchanges the K and K' valleys). 
It implies that 



n T (q,B) = n_ T (-q,-B), 



(29) 



and the relations ( 28 1 follow immediately. The relation 
n T (q,B) — n r (—q,B) that holds if the electric barrier 
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Fig. 4. (a)-(b) Contour plot of the total II c (E,q) and the 
valley II V (E, q) pumping responses for the <5-function magnetic 
barrier |3| in graphene with a gapped spectrum (A = 4El) as a 
function of the canonical momentum q and the Fermi energy E. 
(c) The total Q C (E) and the valley Q V (E) pumped charges as 
a function of the Fermi energy E for different magnetic field 
strengths. 



2 3 4 5 

(E-EJ/E L 

Fig. 5. (a)-(b) Contour plot of the total II c (E,q) and the 
valley II v (E,q) pumping responses for the square magnetic 
barrier (301 in graphene with a gapped spectrum (A = AEl) 



as a function of the canonical momentum q and the Fermi en- 
ergy E. (c) The total Q C (E) and the valley Q V (E) pumped 
charges as a function of the Fermi energy E for different mag- 
netic field strengths. For clarity, the curves corresponding to 
B = 4Bq are multiplied by a factor of 20. 



is left-right symmetric (Uiq = U r o) is not altered by the 
mass term and remains valid. 

In Fig.[4]the pumping responses and the pumped charges 
are shown. The valley pumping response is appreciable at 
low Fermi energy and tends to zero at high energy. It is 
also gradually destroyed by higher magnetic field (note 
that the minimum Fermi energy E m ~ B for high value 
of B). By changing the value of the magnetic field and 
the band gap a high polarization of the valley current can 
be obtained — a regime when current is mainly carried by 
electrons from one valley, so that |Q c (i?)| ps \Q v (E)\. 

We now discuss a square magnetic barrier (with the 
same area) whose position coincides with that of the elec- 
tric barriers (Fig.[TJ;): 



(30) 



In order to apply the transfer-matrix method, we need 
eigenstates in the region with homogeneous electric field 
U(x) — U a and magnetic field B{x) = B. The eigenstates 
can be expressed via parabolic cylinder functions D p (z) 
(sometimes also termed Weber functions) |4*5] : 




To obtain the transmission T T (E, q) and the pumping 
response n T (E,q) we resort to numerical evaluation of 
the transfer matrix. Similar to a (5-function barrier, for 
a device without electric barriers the transmission turns 
out to be the same for both valleys due to the symmetry 
reasons. The valley dependence appears in the structure 
with non- vanishing electrostatic barriers [151116) . 

For a finite width magnetic barrier the pumping re- 
sponse and the pumped charge show similar oscillatory 
behavior, as can be inferred from Fig. [5] An interesting 
feature can be observed in Fig. [SJ;. The total pumped 
charge Q C (E) is a sign-changing function of the Fermi en- 
ergy E, so that at certain values of E, the total pumped 
charge vanishes, whereas the valley charge remains finite. 
In this regime the currents carried by electrons from the 
valleys K and K' flow in the opposite directions: Qk(E) = 
—Qk'(E). They compensate each other (no net charge 
transport), and a pure valley current is generated. It is 
similar in spirit to a pure spin current generation in quan- 
tum pumps. 



ip a (x) 



D p -i(±z) 
±iV2e- 1 D p (±z) 



where 



V2 



£b 
£_b 

L 



£b E — U a + At 



L 



E, 



(E-U a 



A 2 



2EI 



4 Double delta-function and triple square 
magnetic barriers 

(31) The magnetic field profiles created with ferromagnetic strips 
with magnetization perpendicular or parallel to the graphene 
plane are characterized with zero average magnetic field 
(B z ) = [38] . We consider two profiles having this prop- 
erty: a double ^-function barrier with the magnetic field 
in the opposite directions and a triple square barrier. 

(32) Suppose the left gate is non-magnetic and the right one 
has magnetization parallel to the graphene plane. Then, 
the resulting magnetic field can be approximated with 
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the transmission. The vector potential is now left-right 
symmetric (A(x) — A{— xj) rather than antisymmetric, 
and the symmetry of the distributions T T (q) and II T {q) 
with respect to the interchanging of the sign of the canon- 
ical momentum, q —¥ —q, is broken. 

As can be seen in Fig. [6^, the pumped charge Q(E) is a 
peaked function of the Fermi energy (peaks become more 
pronounced as the magnetic field increases). The major 
contribution to Q(E) comes from the quasi-bound states 
between the two 5-function barriers with small longitu- 
dinal momentum k = ki = k r in the leads. Hence, the 
positions of the peaks approximately coincide with the 
energies E n of such states: 



-(e) ! I 


B=BJA ' 

1 B B o 

1 B=4B - 

\ 

A 







E n — E m — El 



2 2 9 
71 71 



n = l,2,. 



(35) 



t 1 2 t 3 4 5 

(E-EJIE L 

Fig. 6. (a)-(d) Contour plot of the pumping response IJ(E, q) 
and the transmission T(E, q) for the double 5-function mag- 
netic barrier ( 33 1 in graphene with gapless spectrum as a func- 



tion of the canonical momentum q and the Fermi energy E. 
(e) The pumped charge Q(E) for a wide ribbon as a function 
of the Fermi energy E for different magnetic field strengths. 



B[x) = 2BL [8{x) - 6(x - L)] z (Pig.jTJi). The correspond- 
ing vector potential takes the form 



A(x) 



(0, 
\2BL, 



x < or x > L 
< x < L. 



(33) 



To model the magnetic field produced by the right gate 
with perpendicular magnetization, we employ a triple square 
profile (Fig. [l^): 



A(x) 



0, 

-B(x + L), 
2B(x-L/2), 
B(2L-x), 



x < —L or x > 2L, 
-L < x < 0, 
< x < L, 
L < x < 2L. 



(34) 



First, we consider a gapless graphene. The structure 
now has a non- vanishing transparency for all energies, but 
when E < E m with a minimum energy E m given by (13), 



which are determined by the simple condition k\L — ixn 
with k\ the longitudinal momentum between the magnetic 
barriers. These values are marked with the black arrows 
in Fig. [6] 

In considering the pumping of the valley current in a 
gapped graphene, we first note that if the vector poten- 
tial is left-right symmetric (A(x) = A(—x)), the Hamilto- 
nian §20§ satisfies the symmetry OH (B)O^ 1 = H(B, Ui <h> 
U r ) with O = TyOyRy. Then, as it was demonstrated in 
Rcf. [TB], the transmission is the same for both valleys, 
Tt(q) = T- T {q). It can be shown that the pumping re- 
sponse Q is also valley-independent, II T (q) = II^ T (q). 
This symmetry is realized, e. g., for a symmetric double 8- 
function barrier B z = BL [8(x + L) — 8(x — L)\ or a dou- 
ble square barrier B z = B [sgn(x + L) — 2 sgn(a;) + sgn(a; — 
No valley current will be pumped in this case. If only one 
of the electrodes has finite magnetization, this symmetry 
is broken. 

The results for a gapped graphene with finite offset of 
electric barriers heights Uio = U r o = — 0.3El are collected 
in Fig. [7j The pumped charges Q C (E) and Q V (E) exhibit 
the oscillatory sign-changing behavior, and similar to the 
structure with a single square barrier considered in the 
previous section, at some values of the Fermi energy a 
pure valley current is generated. 

The dashed lines in Fig. |6^-d and Fig. [7^,-b indicate 
the boundary between propagating and evanescent modes 
under the right electric barrier (in the region < x < 
L). We see that the contribution of evanescent modes 
to pumping is small, but finite. In the case of a gapped 
graphene only these modes participate in pumping at en- 
ergies A < E < E m . 



L)}. 



all modes under the right electric barrier become evanes- 
cent. The analytical expression for the pumping response 
is rather complicated, and below we present the numerical 
results. 

The pumping response, the transmission and the pumped 
charge are shown in Fig. [6] In contrast to a single 8- 
function barrier (cf. Fig. [2] and Fig. [3}, there appear the 
pronounced resonances in the transmission due to the pres- 
ence of a additional magnetic barrier at x — L [10 . The 
similar behavior is exhibited by the pumping response, 
and its resonances are accompanied by the resonances of 



5 Conclusion 

In this paper we have applied the scattering matrix ap- 
proach to the adiabatic quantum pumping in a graphene 
ribbon with a magnetic barrier. By using a 5-function ap- 
proximation for a magnetic barrier profile, an analytical 
solution is derived. We find that a finite magnetic field 
breaks the perfect Klein tunneling so that all propagating 
modes become sensitive to pumping. At the same time a 
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Fig. 7. (a)-(d) Contour plot of the total II c (E,q) and the 
valley II V [E, q) pumping responses for (a)-(b) the double 5- 
function magnetic barrier ( |33[ ) and (c)-(d) the triple square 
barrier (34 1 in graphene with gapped spectrum (A — 4El, 



B = Bo) as a function of the canonical momentum q and the 
Fermi energy E. (e) The total Q C {E) and the valley Q V (E) 
pumped charges for a wide ribbon as a function of the Fermi 
energy E. Constant offset of electric barriers heights is Uio = 
U r0 = -0.3E L . 



magnetic barrier decreases the overall efficiency of a quan- 
tum pump. 

The joint use of a magnetic barrier and band gap 
engineering in graphene gives a way to generate valley- 
polarized currents in graphene-based quantum pumps. The 
parameters of a device can be adjusted such that a pure 
valley current is produced. The pumping is sensitive to 
heights Uio and U r o of electric barriers, and experimen- 
tally it may be easier to vary these parameters rather than 
a band gap value or a magnetic field strength. 

A 6- function and a square magnetic barrier profiles, 
which are employed in this paper and widely used for 
studying the electron transport through magnetic barri- 
ers, can be viewed as simplified approximations to those 
created experimentally by ferromagnetic strips. More re- 
alistic smooth profiles can be analyzed in the same frame- 
work. It is interesting to consider pumping outside the 
bilinear regime where higher valley currents are expected 
to be obtainable. The motion of an electron in a magnetic 
field is also affected by its spin, which should be included 
into the model. These aspects will be analyzed in the fu- 
ture work. 

The considered pump effect might be found useful in 
the field of graphene valleytronics, e.g., as a source of 
valley-polarized and pure valley currents. Experimentally 



they could be detected using, for example, the valley Hall 
effect QI]. 
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